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TRASANA  MEMORANDUM  2-77 


THEORY  OF  STOCHASTIC  DUELS 
MISCELLANEOUS  RESULTS 


1.0  PURPOSE 


This  memorandum  documents  mathematical  derivations  that  have  application 
in  the  use  and  further  development  of  the  theory  of  stochastic  duels. 


2.0  RESULTS 


2.1  Introduction 


Presented  below  are  particular  applications  of  various  aspects  of  the 
theory  of  stochastic  duels  that  the' author  has  accumulated  during  the 
past  20  years.  Some  of  the  material  dates  back  to  his  employment  with 
what  was  formerly  the  US  Army  Operations  Research  Office  that  operated 
through  the  Johns  Hopkins  University  and  later  with  the  System  Develop- 
ment Corporation,  Santa  Monica,  CA. 


2.2  Marksman  Problem  with  Erlang  n Firing  Time  Distribution 
a.  The  Gamma  probability  density  function  is  expressed  by 


f(t)  = t > 0 a>  b > o 

ab  r (b) 


= 0 , elsewhere 


If  the  parameter  b is  restricted  to  the  set  of  positive  integers  b = n, 
(n  = 1,  2,  ...),  the  function  is  normally  called  the  Erlang  distribution 
and  is  expressed  by 


f(t)  = 


tn_1  e~t/a 


an  (n-1)! 


, t > 0 ; n = 1 , 2 , ... 
a > 0 


= 0 , elsewhere 


Tjne&i.S'a, . , aEJifla,  lasfcfldMt  ■ 


I 

% 


This  may  be  evaluated  by  integrating  around  the  contour  of  Figure  1 by 
using  residue  theory  since  the  integral  on  C vanishes  as  R ->  ». 


Figure  1 Contour  of  Integration 

The  jth  residue  is  obtained  by  removing  the  jth  singularity  and  replacing 
u by  u..  Finally,  the  residues  are  surnned  and  multiplied  by  -2tt1  since 

J 

the  contour  is  covered  clockwise.  This  process  provides  the  equation 


-t/a  ^1  t/a  q '/n  el2"j/n 

h(t)  = -4  ; 

aq  n Z— 4 p £ei2nj/n  _ ei2uk/n) 


(7) 


k=o 

m 


where  the  product  in  the  denominator  is’  equal  to  1 for  n = 1. 

c.  If  it  is  now  noted  that  the  mean  of  the  PDF  is  1/na  and,  follow- 
ing usual  convention,  that  the  rate  of  fire  Cr)  is  denoted  by  the  recip- 
rocal of  the  mean.  Equation  (7)  becomes 


np.-e-nrt  V e"rK| 
h(t)  ' Jtl  2-J  "-1  (ei2sj 

<i » *nr  n 


nrtq1/n  e12,j/n 


/0>t.ij/n  i 2rk/n  % 

j.o  n ^ e 1 

k=o 

kft 

With  more  manipulation,  Equation  (8)  may  also  be  written 


n = 2,3...  (8) 


h(t)  - 


npre 


-nrt 


(2<,V")n-' 


n-1  1 /rv 

ei2lTj/n  enrtq1/ 

2-j  "-1 

j=o  H sin  »(k-j)/n 

k=o 
Wj 


» n = 2,3,...  (g) 


EXAMPLE  1:  Letting  n = 1 (remembering  that  the  product  in  the  denomi- 
nator is  1 for  n = 1),  then  * 


h(t)  . p r e'P''4 

and  the  distributicr.  function  (DF)  is 

t 


H(t)  - / h (0  dC  - 1 - e~prt  . 


(10) 


(ID 


EXh.^LE  2.  Letting  n 2 
b(t)  - & 


-2rtf\2r/qt  -2r7q  tl  - -2rt 
— £ p = — sinh  2r  /q  t 


and  the  DF  by  Integration  is 


(12) 


-2rt  r _ 1 

H(t)  = 1 - j^sinh  2r  /q  t + /q  cosh  2r  /q  tj  . O3) 

d.  In  passing,  it  is  also  noted  that  the  DF  of  the  time-to-hit  may 
be  given  in  terms  of  characteristic  functions  as 


f (»)  0 r..e~?l“ldu 

[1  - q $(u)]  u 


04) 


t 

by  noting  that  H (t)  = J h(c)  d£  and  applying  Equation  (51)  of  para  2.7.3  a 

o 

to  this  expression  and  to  Equation  ^2).  The  plots  of  Equations  (11)  and 
(13)  are  shown  in  Figure  2 for  comparison. 

2.3  Tactical  Equity  Duel  with  Erlang  2 Firing  Times2 

a.  For  half  of  the  time,  A sights  B first  and  fires  a round  that 
either  kills  or  alerts  B,  in  which  case  the  duel  proceeds  as  a fundamental 


I 


’ ‘V 


2Ibidi  p 811 


H(t)  *=  1 - [expf-Wr}]/^  [sfn  hZr/q  t cos  hZr/q  c 

Figure  2 Plots  of  Equations  (11)  and  (13) 


I P.'Wll -1-J  T J J ■ Hjpwyip  L I..  Ml rflff/*-*** 


duel.  The  other  half  of  the  time  starts  with  B firing  first.  Thus,3 

PM  » | tPA  + qA  PMf]  + \ qD  PMf  Os) 


where  the  subscripts  refer  to  the  parameters  of  A and  B,  and  P[A]^  is  the 
probability  of  A winning  the  fundamental  duel. 

b.  If  the  known  result  of  the  fundamental  duel  with  Erlang  2 firing 
times4  is  substituted  into  Equation  (15),  the  expression  becomes 


<PArA  - PBrB)t2rA  - <rA  + + 4rArB(rA  4 rB>»A  + (rB  ' rA^ 


<PArA  - PBrB)2  + 4rArB(rA  + rB^ArA  + PBrB> 


(16) 


2.4  Different  Tactical  Equity  Duel 

a.  The  duel  described  in  paragraph  2.3  above  is  modified  for  pro- 
viding a different  approach.  As  before,  each  contestant  fires  one  round 
first  half  of  the  time.  However,  in  this  case,  the  opponent  has  a loaded 
weapon  and  immediately  returns  the  fire  with  one  round,  thus  precipitating 
the  duel  if  both  survive  the  opening  engagement.  From  this  description 


P[A]  = 2 {PA  + qAqB  P^f  } + 2 {qBpA  + qBqA  P^f} 

= PA0  + qB^  + qAqB  P^f  0?) 


v/here,  as  before,  P[A]f  is  the  probability  that  A wins  the  fundamental 
duel . 


b.  If  the  case  of  both  sides  having  negative  exponential  firing 
times  is  substituted,  the  expression  is 


P[A]  = \ PA(1  + qB) 


qAqBpArA 
pArA  + pBrB 


08) 


■J 


Ihid3  Equation  (28) 
Ibidj  Equation  (24) 


6 


'•  '•tV’Wi-’.'S^Wi^i. 


2.5  Duel  with  Random  Initial  Surprise  and  Laplace  Sighting  Times 
a.  The  general  solution  of  this  duel  is5 


P[A]  = 2 h"  f V"u)  Vu^  — 


(19) 


or 


P[A]  = 1 + 2^  /V-o  Vu^  ^y. 


(20) 


where  ./*and  ./are  the  usual  indented  contours  in  the  lower  and  upper  half 

of  the  complex  plane,  and  $^(u),  $B(u),  and  e(u)  are  the  characteristic 

functions  of  the  time  for  A to  hit,  time  for  B to  hit,  and  sighting  time, 
respectively.  A positive  sighting  time  is  an  advantage  for  A,  whereas  a 
negative  sighting  time  is  an  advantage  for  B. 


b.  If  negative  exponential  firing  times  are  assumed  for  A and  B,  and 
if  the  sighting  time  is  allowed  to  be  Laplace-distributed  [g(t)],  then 


, , , _ PArA 

“ p^r^  + iu 


(21) 


*(«>  - ^ 

B PBrB 


iu 


(22) 


g(t)  = 


2c 


Jt-dl 

e c 


- OT  < t < +» 

c > 0 

-»  < d < +» 


(23) 


and 


e(u)  = 


Jdu 


1 + c2u2 


(24) 


5 Ibid , p.  813  et  seq 


ZSM2sjmSiaj**SBL£MaluS*t& 


c.  A typical  PDF  for  the  Laplace  sighting  time  with  d = 0 is  shown 
in  Figure  3.  Note  that  the  mean  is  d,  and  the  standard  deviation  is  /Z  c 


There  are  2 simple  poles  in  the  upper  half  plane,  at  u = ip^  and  u = i 
By  applying  the  residue  theorem,  the  result  is 

p„r„  e-PArAd  P,r.  pBr„  e*'0 


(P(SrB  + PArAJ  0 ' c PM’  2(PBrB. + 71  (pAr 


5(0  = ^ 


c ■ 1;  d » 0 


Figure  3 Laplace  Sighting  Time  Probability  Density  Function 

The  two  cases  that  must  be  distinguished  are  for  positive  d and  negative  d. 
For  positive  d.  Equation  (20)  and  the  upper  half  of  tne  complex  plane  are 
used.  Thus, 


_ , ! r PBrB  e 

P[A]  = 1 + Jpr  / 

2 J (pArA  + iu^pBrB  “ iu)0  + 


du 

c2u2)  u 


j £>.■ A j?-,  k i,<v. 


d.  Two  further  simplifications  may  be  noted.  If  it  is  supposed 

that  A always  has  the  sighting  advantage  and  that  g(t)  = 1/c  e-t^c, 
t > 0,  c > 0,  and  g(t)  = 0 elsewhwere,  and  e(u)  = 1/(1 -icu),  it  is  easy 
to  calculate  as  above  and  show  that 


pm  . + PArA  PBrB...  _ (30) 

lPArA  + PBrB)(1  ' P8rBc)  (PArA  + C)(pBrB  ' c1 


ffr 

Similarly,  if  B always  has  the  sighting  advantage,  with  g(t)  = 1/c  e , 
t < 0,  c > 0,  and  g(t)  = 0 elsewhere,  and  o(u)  = 1/(1  + icu);  then 


PC  A]  = 


PArA 


PaVW(1  + p«rRc)  ‘ 


In  this  case,  obviously,  P[A]  has  been  reduced,  compared  to  the  funda- 
mental duel,  by  a factor  of 


pBrBc 


2.6  Fundamental  Duel  with  Erlang  Firing  Times 


a.  In  this  application,  the  fundamental  duel  is  examined  where  A 
has  an  Erlang  n firing  time  PDF,  and  B has  an  Erlang  m PDF.  This  means 
that 


tn-l  -t/a 

f„(t)  = K :-e - 


an  (n  - 1)! 
= 0,  elsewhere 


t > 0,  a > 0 
n = 1,  2,  ... 


with  4»a(u)  = - — -=•  » 

A (1  - iau)n 


.m-1  -t/b 

fR(0  " ^r~~7~  • t > o,  b > 


bm  (m  -1)!  m * 1,  2,  .. 
= 0,  elsewhere 


with  $B(u)  = - 


(1  - ibu)r 


rr*\TtXZ,+‘ v W-y" ih'y^ifj. a s**  *•%  j*3 


From  previous  work,  we  have6 


PM  = J-  f — 

2*’  J 10  + lau)"  - qj 


du 


qA]  [0  - ibu)ra  - qB]  u 


~ (34) 


which  may  be  rewritten  as 


P[A]  = 


PaPb 


l-itD 


r 1 d 

jp  + i^u)"  - qjj'ltoi)"-  (£)\]“ 


du 


2ni 


(35) 


There  are  m simple  poles  in  the  lower  half  plane,  at  \ - - r(1  - qVm  el2,lk/m)) 


uk  " ' p 1 ‘ qB 

k = 0,  1,  ...,  m-1.  Hence,  Equation  (35)  may  again  be  rewritten  as 


if 


Dr.,  pApB  f 1 V du 

1 1 " M J [(1  * iau)n  - q„)  ‘ “V  [„  + i{l  - qf"  ei2=fc/n)J  » ‘ L ‘ 


k~o 


• th 


b.  The  j singularity  is  now  removed,  and  u-  is  substituted  for 

J 


♦ th 


all  u’s  to  get  the  j residue.  The  residue  theorem  is  then  used, 
multiplying  by  -Zui  and  summing  over  all  residues.  Finally,  it  is  noted 
that  r^  = 1/na  and  r^  = 1/mb.  Thus 

m-l 


PM 


.!£»  V ~T - 

■f"’"  % {>  - 4'° 


(37) 


WJ 


where  it  is  noted  that  the  product  term  in  the  denominator  is  1 if 
m = 1.  One  can  immediately  check  this  result  by  substituting  m = n = l7 
and  by  letting  m = n = 2. 8 Equation  (37)  above  may  be  written  in  the 
following  form  which  may  be  somewhat  simpler  to  evaluate. 


PM 


- PAPB  V 

(2qB/0)K'  {[1 


i . n5 


q,(l  - qj/n  e12,JA=)}n  _ qJ(1  . qVo  'iZ.j/C)  o-l  j1n  j(k  . j) 

• rA  t . « 

VI  i 

and  again  the  product  term  in  the  denominator  is  1 if  m = 1. 


(38) 


Ibidj  Equation a (14)  and  (20) 

7 


Ibidj  Equation  (0) 


Ibid , Equatioti  (24) 


12 


jaeskaas 


i.  jd  aidfla  ■ aaeflBajBa^  aca 


- .<  -*  driAhcj* aca  ■■*  saSaaesJI 


■ - -*-\  "*r-\  ***»<o.*'  ,<i 


• r*  •» *»  . ♦«•  • 

J Giut  fA(t)  fg(t)  fc(t)  dt  **  ^2  J ♦*(«-«)  f J ❖gCw-v)  «c(v)  dvj  dw  (43) 

- <o  - CO  \.w  / 


This  result  comes  about  by  applying  Equation  (42)  to  the  left-hand  side 
twice. 


N.B.  These  theorems  do  not.  require  the  functions  f(t)  to  be  pdf's.  They 
are  general  for  any  integrable  functions. 


2.7. 2 Properties  of  Characteristic  Functions  of  Positive  Random  Variables 
a.  4(0)  « 1 (44) 

PROOF:  From  Equations  (39)  and  (40) 


00 

•/ 


uo)  t 


4(0)  = / e""1'  f(t)  dt  - J f(t)  dt  = 1 

o o 


00 

■/■ 


b. 


1 <»(u)  | <1-,  (Imaginary  u > 0) 


(45) 


This  implies  no  singularities  in  the  upper  half  of  the  complex  plane. 
PROOF:  Let  u = v t iw,  then 


CO 


Mu)  » J eivt  e"°*  f(t)  dt 


or 


CO 


Mu)  | 5 J e"1*  f(t)  dt  , |e£ut  | <1 


and 


CO 


4(u)j  - J f(t)  dt  = 1 

o 


e ot  I i 1 


■p 

§ 


1 


3 

1 


I 


i 


<5 

l 


l 


3$ 


IS 

I 


i 

w 


i 


$2 

|2 


1 


<3 


& 


■Q 

% 

1 


J 


to 


% 


for  w > 0 , i.e.,  positive  imaginary  u. 


>■  mmmatmMtaiiimmsm 


|<K-a/|  < 1 (Imaginary  « < 0) 


This  implies  no  singularities  in  the  lower  half  of  the  complex  plane. 


PROOF:  Same  as  for  subparagraph  b above,  with  negative  u where  final  step 
depends  on  jewtj  < 1 for  negative  w. 


f(t)  a differentiable 
function  of  bounded 
variation.  Imaginary  u > 0, 
k = positive  constant. 


This  implies  that  <f>{uj  diminishes  as  1/R  in  the  upper  half-plane  where  R 
is  the  radius  of  a semicircular  path  of  integration  in  the  complex  plane. 


PROOF:  Integrate  Equation  (40)  by  parts  to  obtain 


CO 

- i.mi  + £ /V 

u u J 


f'(t)  dt 


This  makes  use  of  the  fact  that,  necessarily,  f( 0 = 0.  Now,  for  u complex 


and  t > 0,  JeTU  | < 1 in  the  upper  half-plane.  Thus, 


CO 

l*wi  ^ + w/lf,(t)l  dt 


Now,  assuming  a differentiable  function  f(t)  of  bounded  variation,  then 


CO 

J Jf' ( t) j dt  is  bounded;  therefore,  \§(u)\  5 where  k is  a posits 


constant  and  imaginary  u is  positive. 


f(t)  a differentiable 
function  of  bounded 
variation.  Imaginary  u < 0, 
k s positive  constant. 


This  implies  that  <t>C-uJ  diminishes  as  1/R  in  the  lower  half-plane. 


-•i-  wk+i 


I ! 


^^^r<Tr-jr”’^'~"T"?r.^-; 


PROOF:  Same  as  for  subparagraph  d above  using  negative  u.  In  this  case 
for  u complex  and  t > 0,  je"iufc  j < 1 in  the  lower  half-plane. 


f. 


\$(u-w) | < 1 (Imaginary  u > 0) 

(Imaginary  w < 0) 


(49) 


This  implies  no  singularities  in  the  upper  half  of  the  complex  u plane 
and  the  the  lower  half  of  the  complex  w plane. 


PROOF:  Let  u = C + £n  and  w = x + iy  so  that 


co 


$(u~w) 


- /•**' 


■”ew’  f(t)  dt 


and 


CO 


k! 


u-u) I < f 


e<V-o>t|f(t)  dt 


co 


■-I 


f(t)  = 1 


if  y- n < 0 for  all  y, n.  This  implies  y is  negative  and  n is  positive  or 
that  Imaginary  u positive  and  Imaginary  y is  negative. 


| u-u  j 


f(t)  a differentiable 
function  of  bounded 
variation.  Imaginary  u > 0, 

Imaginary  w < 0,  and 

k = positive  constant. 


(50) 


This  implies  that  $ (u-u)  diminishes  as  1/R  in  the  upper  half  of  the 
u plane  and  ‘in  the  lower  half  of  the  w plane. 
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PROOF:  Integrating  the  integral  form  of  $(u-w)  by  parts,  again  using 

f(co)  = 0 


= £iI°U  ±-  / eiMtf(t]  dt 
o 


Now,  |eir“‘u;tJ  < 1 for  Imaginary  u positive  and  Imaginary  w negative. 
Therefore , 


i i^r+jH-W  /|f,(t)!  dt 


- and  again,  if  f(t)  is  a differentiable  function  of  bounded  variation 


\$(v.-w)\  < where  k = positive  constant. 

“ |«-wj 

2.7.3  Additional  Theorems 

a.  a 

Jeiutf(t)  dt  = ^ J * M>(1  e-^-  dw  (51) 

O - CO 


t 

i 


Note  that  for  w = 0,  an  expression  is  also  given  for  the  distribution 
function  of  a random  variable  in  terms  of  characteristic  functions. 

PROOF:  Using  Equation  (40) 


-zk) 

— co 


..'-*---1  r-L  d cbm.  »-  - j>a<-  --a*'*. 


The  desired  result  is  now  obtained  by  replacing  w with  (w  + u). 


♦ co 


/ e<UC  f(t)  dt  = *(u>  - Hf  / ^f-  - e"“)-  * (52) 


- CO 


Mote  that  for  u - 0,  this  provides  an  expression  for  the  complementary 
distribution  function  in  terms  of  characteristic  functions. 


PROOF: 


CO 

f> 


r 41 

e'“cf(t)  dt  - J eiut  f(t)  dt  - J eiut  f(t)  dt 

a O o 


= »(U)  - J,  T. 

Zin.  J w 


dw 


where  we  have  used  Equations  (40)  and  (51). 


c. 


/V(/f{?)  dE)dt 

o t 7 ♦ 


(53) 


This  is  the  characteristic  function  of  the  complementary  distribution 
function. 

PROOF:  Integrating  by  parts  results  in 


to  « 


fe6*  ( [ f(t)  dE  1 dt  - - 4-  + , ±£hJ_-_1 

J \ J J 1,U  1U  %U 

A ’ » • 
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d. 


r-  » 


«v  v ~|  4 «0 

/fAtt)  /vodcldt-ijv^ii 

o L t J -o 


r*»r«j  - 1] 


du 


(54) 


PROOF : Pre vi ous 1 y gi ven . 1 1 


dj  1 du 


(55) 


f. 

T w 


+ 05 


+ «S> 


M/v.i  *)  * ■ ijK-Kr" ' ‘y-1  - ” *) 


<&))  <2w  (56) 


PROOF:  Previously  given.13 

g. 


ffA(t>  (fVc>  d(j  dt  “ 


r t.  _«a 


£AfK  - w;[<Bfu;  - i] 


w 


&>)<?«  (57) 


n 


Trevor  Killians  cold  Clinton  J.  Ancker , Jr,  "Stochastic  Duels, " Operations 


Research , Volume  11,  No.  5,  October  1963 , pp  603-817,  Equation  (18). 
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C.  J.  Anckcr,  Jr,  "Stochastic  Duels  of  Lir.tiied  Time  Deration'," 


CORS  Journal,  Volume  4,  No.  2,  July  19CC,  pp  69-61,  Equation  (S). 
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Ibid,  Equation  (28) 
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Now,  let  - p to  provide 


= _L 

2-ai 


+ CO  ^ 


and  by  Equation  (40) 


, tul-auK^M  - 1] 


T co 

/'*«*>( 


at  + b 


J e~iOT  - 1 ^ j\(U  ~ atf)  ~ 1]  ^ ^ 


du  (61) 


PROOF:  Let  t = ^ , then 
a 


at  + b 


/fAtt)  ( / fBw  d0  dt)  * 

n ^ . n 


at  + b 


and  replacing  the  inner  integral 


at  + b 


- [ t ( f i-  [ 
-J  fA  V a;\  2tt  J 


1 C*B(«)  - 13 


’)r 


at  ♦ b 


1 f ^2  ^ /"f  /n-A  c-wodn\ 

'*  2ttt-  J v \ / ^ \ a J a * 


PL . * •,& y. fcfc* I'&rvStr-fC’jr  , *-  . 


Now,  let  = p 


and  using  Equation  (51) 

. 1 " 13  -a*  / 1 f *A(“  “ GZj)tl  " e~lUXl 


t i r 9p\u  - com 

(Zrf  J u 


d«  } cb 


i - u w fV<“'  - 1]  y*  - ^ ^ 

"4 i*- J 5 e II £ <!u)da 

_*rrt  — 00  * 


which  provides  the  form  of  Equation  (61)  if  the  order  of  integration  is 
reversed. 


w co 

/v«(  /yodel** 

T it  ♦ b 


f“  e-**  ♦.(-»)[*„(«)  -1] 


= zli/ 


1 f [e'iut  - 11  / f $pSu  ~ ^ e rW&C4>B(u)  - 1]  \ 

•drj  *H"  ](J~ 5— J 5 d“)d“  («) 


PROOF:  The  integral  may  be  written  as 


to  co  M 

JfAlt}(  / fB<  0 d«  ) « - / fA(t)  ( f fBU)  dc  ) dt 
0 ° «U 


= Equation  (60)  - Equation  (61) 


i._/>fA«0  (/ fB(e>  is)dt  . ^ f A (63) 

0 ' t w 

PROOF:  First  apply  Euation  (42)  and  then  Equation  (53)  to  the  left-hand 
side  of  Equation  (63)  to  obtain'  the  right-hand  side  immediately. 


(64) 


PROOF:  First  apply  Equation  (55)  and  then  Equation  (63)  to  the  left-hand 
side  of  Equation  (64);  with  slight  rearrangement,  one  obtains  the  right- 
hand  side  immediately. 


i 7 - 1] 

’ 4.W  4.« 

r [^(u-u)  - 13/  f ^{-^[^(w-y)  - 1] 

8 J v 

***W 

J ( i>-v ) \ J (u- v) 

— W -M  ' 

cfo 


dy 


(65) 


PROOF:  Apply  in  succession  Equations  (42),  (43),  and  (53)  to  the  left- 
hand  side  of  Equation  (65),  and  with  some  rearrangement  of  terms,  the 
right-hand  side  is  readily  obtained. 
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'j'u  - •njp'v® 
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Wj 

I 


I 

I 

l 

as 


The  substitution  now  made  is  -u  = x + -Ly  which  is  an  analytic  continuation 
of  the  real  variable  u into  the  complex  plane.  This  substitution  provides 


P[A]  = |7  J fit)  dtj 


txt  -srt 
e e 

x + iy 


(dy-idx) 


f(g  sin  xt  - y cos  xt)dx  + (x  cos  =t  + v sin  xt)d^~] 
s2  + y2 


~yt 


Ux  cos  xt  + v sin  xt)dx 


* IU  cos 


X1  + y2 


(68) 


c.  Consider  the  two  integrals  within  the  brackets  in  Equation  (68); 
specifically: 

(1)  By  integrating  along  any  line  parallel  to  the  x-axis  at  some  y, 
e .g. , y (which  may  be  +,  -,  or  zero,  thus  including  the  x-axis),  dy  = 0, 

and  the  integrand  of  the  real  part  is  even  and  the  integrand  of  the  imagi- 
nary part  is  odd. 

(2)  By  integrating  along  lines  parallel  to  the  y-axis  at  some  fixed  x, 
e.g.,  xq,  dx  = 0,  and  the  integrands  are  neither  even  or  odd.  However,  if 

two  lines  are  considered,  one  at  +x  and  one  at  -x  , it  is  noted  that  the 

integrand  of  the  real  part  along  -x  is  the  negative  of  the  integrand  of 

the  real  part  along  +x  , i.e.,  the  integrand  is  odd  relative  to  the  x variable, 

and  the  integrand  of  the  imaginary  part  is  exactly  the  same  at  +xo  and  -xQi 

i.e.,  the  integrand  is  even  relative  to  the  x variable. 


% 

I 


d.  The  result  of  the  above  is  that  P[A]  may  be  evaluated  along  the 
contour  shown  in  Figure  8.  The  property  described  in  c(l)  above  means 
that  the  imaginary  part  along  Cj  is  equal  to  the  imaginary  part  along  CG 

except  that  the  signs  are  opposite;  consequently,  the  sum  is  zero  while 
the  real  parts  are  equal.  The  same  applies  for  the  integrals  along  C3  and  C4. 


<-xo.0)  (0.0)  (+xQi0)  Cc 


CV^o* 


Cs  (o.-y0)  ) c*. 


(4x0,-y0) 


Figure  8 Numerical  Integration  Contour 


e.  From  the  property  described  in  c(2)  above,  the  imaginary  part 
of  C2  is  again  canceled  by  the  imaginary  part  of  C5>  and  the  real  parts 

are  equal.  Thus,  the  only  thing  that  must  be  done  is  to  integrate  the 
real  part  of  the  integrand  along  C4,  C5,  and  Cg  and  multiply  the  sum 

by  2 to  get  P[A].  This  can  significantly  reduce  the  amount  of  numerical 
integration. 

f.  One  should  also  note  that  exactly  the  same  argument  holds  for 
Equation  (19)  or  for  any  integrand  that  is  a product  of  characteristic 
functions  v/ith  at  least  one  positive  argument  and  at  least  one  negative 
argument  divided  by  u.  In  addition,  one  may  speculate  that  -y  and  +xQ 

should  be  chosen  so  as  not  to  be  too  close  to  the  singularity  at  the 
origin  or  to  the  nearest  singularity  in  the  lower  half-plane.  This 
would  avoid  steeply  varying  integrands  and  should  reduce  numerical  inte- 
gration difficulties.  Also,  one  may  simply  continue  C4  to  «,  thus  elim- 
inating C5  and  C6. 

2.8  Generalizations  of  Some  Results  by  Thompson 

Thompson  has  derived  some  results  for  incorporating  reliability  into  the 
theory  of  stochastic  duels.14  However,  he  has  not  simplified  them  by  using 
characteristic  functions.  This  has  the  disadvantages  of  leaving  the  solu- 
tions with  one  or  more  infinite  sums  of  weighted  convolutions  of  pdf's  to 
be  integrated.  These  operations  can  only  be  closed  for  very  special  pdf's 
and  thus  are  of  limited  utility.  By  using  characteristic  functions,  these 
sums  can  always  be  closed  to  a simple  form,  and  the  number  of  multiple 
integrations  to  be  performed  is  reduced  by  at  least  one.  Finally,  there 
exists  a characteristic  function  for  every  pdf  and  if  it  is  unknown,  can 


David  Thompson , Development  of  Analytical  Models  of  Battalion  Task 
Force  Activities , Systems  Research  Laboratory , Department  of  Indus- 
trial Engineering , University  of  Michigan,  Ann  Arbor,  Michigan, 
Report  Number  SRL19757  FR70-1  (U),  Part  F,  Chapter  1,  Editors:  Seth 
Bonder  and  Robert  Farrell. 


be  approximated  from  real  data  for  practical  applications.  For  these 
reasons,  we  now  simplify  some  of  Thompson's  results  as  follows: 

a.  In  the  case  where  either  side's  weapons  may  fail  to  function 
after  some  period  of  time  and  no  withdrawals  are  permitted,  Thompson 
provides  the  following  expression  for  the  probability  that  side  A wins:15 


CO  m ^ to 

P[A]  = J hW/, rAU)  d Ah  - / hg m([ drA  dx 
o \ L o x / - 


dt 


(69) 


where  hA(t)  and  hg(t)  are  the  usual  pdf's  of  A's  and  B's  times  to  a kill, 

and  rA(t)  and  rg(t)  are  A's  and  B's  pdf's  of  times  to  weapons  failure. 

All  are  independent  of  each  other.  We  now  rewrite  Equation  (69)  in  a 
more  convenient  form  for  our  purposes,  i.e.. 


+ 


rAU)  dt 


hB(x) 


rg(n)  dn 


dx 


dt 


(70) 


We  now  apply  Equations  (54)  twice  to  the  first  expression  on  the  right- 
hand  side  of  Equation  (70)  and  Equation  (64)  to  the  second  expression  to 
obtain 


P[A] 


, rcAt-*)[?A(to  - 13 

2 r.i  J u 


1 ^ 
1 -2ST  J u 


do 


hfuf 


«B(-u-v?[rB(p)  - H 4t> 


\f  r*«A(u-o)[TA(«)  - 1]  \ 

* VW  » / (7,) 


where  $(•)  is  the  characteristic  function  of  h(t)  and  v(*)  is  that  of  r(t). 
Since  all  of  the  pdf's  are  for  positive  random  variables,  we  may  now  use 


15 


Ibid  p582 } Equation  2 


>S  ■4r»--57  rf  ‘^hMii^st*  ^ _ 


the  properties  of  positive  RVs  given  earlier  to  further  simplify  these 
results  to 


r , 1 f* a<-u)  fA^ 

P[A]  « j- J -A-  --^-  du 


1 f «B^w)  Vu)  , 

1 "-ttj  ^ 

L 


W 


8ir  3f 


i r i / r V-"-^  . \ / r *a («-»)  V“)  \ . . 

i?  — v — -d^y  a...  dJjdK>  (72) 


V L 


The  corresponding  expression  for  a draw,  P[AB],  is  provided  by  Thompson16 
as 


P[AB] 


« ^ O ^ 

= J ^*4  f hAU)  ds  di*  J rB(t)f  f hB(*)  dxl  dt-  (73) 

o L t o *-t  J 


Applying  Equation  (54)  to  each  half  of  this  product,  we  have 


p Cab] 


1 *f  *a(-)I**W  - U TR(-u)[>R(t>)  - l] 

" >*V 


4n2 
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dit 


•/- 
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dw  (74) 


which  is  now  further  simplified  using  the  properties  of  positive  random 
variables  to 


P[AB]  = - 


1_  f YA("u)  *a(u)  dy  . *B(u>) 


47T2 


/■ 


u 


/ 


w 


cb  . (75) 


Thompson  provides  an  example17  where  all  the  pdf's  are  negative  exponentials. 
When  these  pdf's  are  used  in  Equations  (72)  and  (75),  we  check  out  with  his 
results  except  for  an  error  in  his  expression  for  P[A].  He  should  have 
added  a in  the  last  parenthesis  of  the  denominator. 


16 Ibid , p5B2j  Equation ' 3 
17 Ibid , pr>83 
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Again,  we  apply  Equation  (65)  to  each  term  of  the  right-hand  side  of 
Equation  (79);  thus 


P[A3] 


1 f CoB(»)  - 13  r *r  C V**-")  “ *3/  f 

ZnKJ  v J u-y  [J 


lt>Aa-v)  - 1]  / /•r.t-alhnM  - 1] 


Av/l8v 

U-tf 


■) 


du  ) cL) 


di> 


j r CCa(u'v) : 13 / ' - dnUI 

" B«%J  v J v-v  \J  u-u  / j 


- du  1 du  du  (80) 


which  again  reduces  to 


L \ LL 


tb)  dy  (81) 


and  the  indentations  are  the  same  as  for  Equation  (78).  Again,  Thompson 
provides  an  example  using  all  negative  exponential  density  functions.20 
When  these  assumed  functions  are  inserted  in  Equations  (78)  and  (81)  and 
are  integrated,  the  results  check  with  his. 

c.  In  the  next  duel,  we  have  the  more  difficult  case  where  a duelist 
withdraws  when  his  weapon  fails  but  only  after  he  discovers  it,  which  is 
on  the  next  attempted  firing.  For  this  case,21  Thompson  provides 

PEA]  = /;,)(/ rA(c)<J^[^yrB(2/)d^  J *Bf!  ^d*  dt  (82) 

We  note  here  that 


V n 

bRU) 

PB  B 


and  shall  hence  forward  use  that  fact. 


20IbidJ  p586 

2^Ibid}  p588,  Equations  (6)3  (7),  and  (8) 
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Now,  by  writing  out  the  entire  first  term  of  Equation  (83)  and  reversing 
the  order  of  integration,  we  get 

LHS  of  the  first  term  of  Equation  (83)  = 
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PR  En3£ 


j mi- 


c-B(y+u-p)i>B(f>)-l] 


!*)  (/t''t 


(v+u)t  -lut 


)bA(t) 


(87) 


The  inner  integral  of  Equation  (87)  is  from  the  application  of  Equations 
(42)  and  (53),  therefore 


Now,  substitute  Equation  (88)  into  Equation  (87)  and  then  Equations  (87) 
and  (84)  into  Equation  (83),  thus 
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(89) 


Equation  (S3)  is  now  further  reduced  using  the  properties  of  positive 
random  variables  to 


P[A] 
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(90) 

II  L ' 

This  expression  again  checks  Thompson's  example22  using  all  negative  expo- 
nential pdf's.  The  prior  results  were  based  on  a failure  mechanism  that 
was  a random  function  of  time.  In  what  follows,  the  failures  only  occur 
at  the  instants  of  weapon  firings. 

d.  The  no-withdrawal  case23  is  identical  to  the  author's  random 
ammunition-limited  duel,  and  the  characteristic  function  form  of  the  solu- 
tions is  given  in  the  referenced  paper  and  will  not  be  repeated  here. 

e.  Similarly,  in  the  same  paper  mentioned  above,  the  characteristic 
function  form  of  the  solution  to  the  duel  in  which  a contestant  withdraws 
the  instant  a failure  occurs  is  also  given  and  will  not  be  pursued  further. 

f.  The  case  in  which  a contestant  whose  weapon  has  failed  on  the 

n^  round  but  who  only  discovers  that  fact  on  the  (n+l)st  round  and  then 
withdraws  is  provided  by  Thompson24  as 


(91) 


Thompson  assumes  the  probability  of  failure  on  any  given  round  is  fi'.ed 
at  u and,  therefore,  p + q + u - 1.  Again,  noting  that 


£qBfB  (x)d*  = j7hB(.x) 


- qB 


22  Ibid,  p589 
2 3 Ibid,  pS94 
2 ''Ibid,  pp  597 , 598 


we  rewrite  Equation  (91)  as 
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P[A)  =yhA(t)^y'fB(x)dxj  dt  + y*hA(t^y*hB(x)f  fwj 


dr 


dt.  (92) 


Using  Equation  (54)  in  the  first  integral  and  Equation  (61)  in  the  second, 
we  have 
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In  the  second  integral,  we  now  reverse  the  order  of  integration  and  inte- 
grate out  t to  get 
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Now,  taking  advantage  of  the  properties  of  positive  RVs,  Equation  (94) 
becomes 
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The  expression  for  a draw25  is  composed  of  two  terms,  each  of  which  is 
similar  to  P[A]  and  the  development  is  almost  identical;  therefore,  we 
shall  merely  write  down  the  results,  i.e.. 


P[AB3  = ■/.*<->*«**  * foB(-u)^M<iu/u 


pTO^j/  w 1 J u 

l vu 
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dw ) <±x>  « (95) 


Equations  (95)  and  (96)  both  check  with  one  of  Thompson's  examples26  where 
again  he  uses  negative  exponential  pdf's. 

g.  The  remainder  of  Thompson's  paper  deals  exclusively  with  negative 
exponential  pd47’ . from  the  beginning,  and  all  the  results  are  consequently 
so  specialized  o that  function  that  no  further  simplifications  through 
the  use  of  chai  • teristic  functions  are  possible. 


h.  In  a paper  published  in  the  Naval  Research  Logistics  Quarterly, 
Thompson  extends  the  situation  in  paragraph  f above  to  the  case  where  each 
contestant  has  an  independent  probability  distribution  on  the  pos.-ibility 
of  a failure  on  the  nth  round  (and  thus  a withdrawal  on  the  n+lst  round).27 
He  provides  the  following  as  solutions: 


If  a-  = P [A's  weapon  fails  on  round  j+1]  , 

1 

Bj.  = P [B's  weapon  fails  on  round  k+1]  I 


and 


CD  CD 


(97) 


25Ibid,  pSB9 
2^Jbid,  pp599,  699 

27 David  E.  Thompson , Stochastic  Duels  Involving  Reliability,  The  Naval 
Research  Logistics  Quarterly,  Vol  19,  No.  1,  March  1972,  pp  147-148, 
Equations  (6)  through  (11). 
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